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Eppur si muove - Galileo 

Abstract. In this paper we construct a family of exceptional Stewart-Gough 
Mechanisms, whose motion curves are algebraic curves of genus 7. Up to trans- 
lations, rotations and dilatations this family of mechanisms is 13-dimcnsional. 

Introduction 

We call a mechanism exceptional if its degree of freedom exceeds the dimension 
of a general mechanism of the same underlying topological type. In this paper we 
establish the existence of exceptional Stewart-Gough platforms by a combination 
of algebraic geometric and number theoretic methods. 

A Stewart-Gough platform consists of two rigid bodies, one of which is grounded, 
linked to each other by six telescoping legs in ball joints at both ends of each leg 
(see Figure 1). 



Figure 1. A general Stewart-Gough platform, Sommese- Wampler 
[SW05j . pp. 107. 

Usually, for given leg lengths and rigid bodies, there is only a finite number of 
ways to assemble the parts. In robotics, the ungrounded body is moved around by 
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changing the leg lengths. Since the group of movements SE{3) is six-dimensional 
and since, in general, when its length is locked, each leg places one constraint on 
the motion, six legs is the minimum number to rigidly constrain the ungrounded 
body. As the leg lengths are changed, the ungrounded body can be moved through 
a region corresponding to an open part of SE{S). 

In this paper we are interested in a different situation: We wish to find a 
configuration of leg lengths and bodies such that the mechanism moves without 
changing the leg lengths. We call such a mechanism exceptional. 

Apart from rather trivial examples, one family of exceptional mechanisms is 
well-known and thoroughly studied, the Griffis-Duff'y platforms |GD93j , Husty and 
Karger |HK97j . The family of Griffis-Duffy platforms forms a unirational variety, 
Karger |K03j . In this paper we ask whether there are families of exceptional 
Stewart-Gough platforms of different kind. In our main result, we establish the 
existence of a family of mechanisms whose motion curve in SE(3) is a smooth 
algebraic curve of genus 7. In Figure 2, we display the coupler curve, i.e., the curve 
traced out by a coupler point fixed to the ungrounded link under the motion. 



Figure 2. A Stewart-Gough platform with a coupler curve of de- 
gree 12 and genus 7. 

Our methods, we believe, are quite novel in this kind of research. We first 
approach the question by searching for mechanisms defined over a finite field. We 
then utilize the algebraic-geometric properties of the discovered mechanisms, to get 
a construction over some algebraic number fields K, i.e., fields K with Q C if C C 
of degree [if : Q] dimq K < (x. Finally, to establish the existince of mechanisms 
defined over M, we prove that some of these fields K have embeddings if ^ K. 
The final existence result can be verified also by numerical methods as was pointed 
out to us by Charles Wampler. 
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1. The algebraic setup 

A movement in the euclidean 3-space is a map 

with an orthogonal matrix A = [aij) £ 50(3) and a translation vector b G M'^. We 
denote the group of movements by SE{2,). 

A position of a Stewart-Gough platform is determined by the positions of the 
six legs. Let g*^') G M'^ for z = 1, . . . , 6 denote six points on the grounded body, and 
let p*^*) S be the end points on the moving body. Consider the equations 

eg(*) : || p« - f =|| ApW + 6 - g« for t = 1, . . . , 6, 

where || x |p= {x,x) denotes the square of the euclidean norm for a; G M"^. Let 

$ = y(eg«,i= L...,6) C X SE{i) 

be the complete parametrized family of Stewart-Gough mechanisms and their as- 
sembly positions. Consider the natural projections tti : $ ^ A'^^ and 7r2 : <f> 
SEii). 

Then Cq — ■K2{t^i^ {ct)) is the asssembly space of the mechanism specified by 
a — {p^^\q^'^\ . . . ,p^^\q^^^) G A'^^. For generic a, Ca is a finite set of isolated 
points, i.e., dimCc = 0. We are interested in finding exceptional mechanisms a 
such that dimCa >0. In particular, when dimCa = 1, we call the 1-dimensional 
part of Ca the motion curve of mechanism a. For such a the assembly space consists 
of the motion curve and possibly a finite number of isolated points. 

The equations A = E^ and det A = 1 defining SE{3), and the equation eg^'-' 
make sense as polynomial equations over arbitrary fields K. We call a G A^^{K) 
such that the solution space in Ca C SE(S,K) is one-dimensional, an exceptional 
mechanism defined over K. Here K denotes the algebraic closure of K. For fixed 
a G A^^{K) one can determine the dimension of the solution space in SE{3,K) 
with a single Grobner basis computation over K. Macaulay 2 |M2j . which we 
used, and Singular |GPS01| are good computer algebra systems to perform such 
computations. 

The complete computer algebra code used is documented online in [GS09| . 

2. Experimental exploration 

Let M C A'^^ be a codimension c irreducible component of the constructible 
set of exceptional mechanisms. Heuristically, for the finite field Fg with q elements 
the ratio of points on M is approximately 

iA36(F,)r'' ■ 

By the Weil formula, this equation is asymptotically correct for g — > oo in case of 
an absolutely irreducible component M defined over Z. 

We however plan to make experiments over very small finite fields. We im- 
plement a script for parallel search for the computer algebra system Macaulay 2 
IM2] running on a cluster with 24 nodes with 2.4 GHz each node. To compute 
the dimension of the solution space in S'£'(3,Fq) for one platform takes about 0.3 
seconds. Thus, for very small fields k = Fp with p G {3, 5, 7} we are able to test 
p^p platforms, — 12, A^5 = 9 and ~ 7, each in a maximum time of about 3 
hours. 
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We have reasonable hope to find points on Af (F3) for components M C A^^ of codi- 
mension < 12. To give an idea of such exploration we give the result for a typical 
test as described above. We list the dimension, the degree and the total number 
of the discovered exceptional mechanisms of that degree. Here we use the naive 
embedding SE{i) = 50(3) x c A^^ c P^^. In Table 1 below, the degree reverse 
to the degree of Ca with respect to the naive embedding. To get an interpretation 
for the total numbers, we also print the hypothetical codimensions, which are only 
correct, if the family is absolutely irreducible. In case of a family consisting of 
several components, say fc, of the same codimension the values hcodim are more 
likely approximations of the product h — log k and the codimension codim. 



dim 


deg 


F3 


hcodim 


F5 


hcodim 


Ft 


hcodim 


1 


1 


2817 


4.77 


861 


4.80 


60 


4.90 




2 


6223 


4.05 


1553 


4.43 


85 


4.72 




3 


425 


6.49 


20 


7.14 
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2256 


4.97 


542 


5.09 


29 


5.27 
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383 


6.59 


35 


6.79 


2 
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1083 


5.64 


580 


5.05 


66 


4.85 
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49 


8.46 


2 
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236 


7.53 
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8.00 
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12 


9.74 


1 










10 


23 


9.15 
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12 


43 


8.58 
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8.14 








14 


5 














15 


1 














16 


2 
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1 














24 


3 














28 
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32 
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36 






1 
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4 


8 














6 


9 




4 


8.14 








8 


1 












sum 


13484 


3160 


242 



Table 1. Experimental data 



In the following we investigate the examples of degree 6 or 12 further, leaving 
all other cases to future research. The main reason to choose these cases is that 
they have very nice Betti tables, see below. To get projective curves, we work in 
the projective closure SE{3) C P^^ of SE{3) C A^^. 

We first remove zero-dimensional components of the assembly space, i.e., we 
compute the homogeneous ideal of the one-dimensional part of the assembly space. 
This is a standard task in computational algebraic geometry, which for example can 
be performed using primary decomposition. Since this is fairly expensive and more 
than we want, we actually use an alternative approach built around a probabilistic 
algorithm, see our documentation |GS09| . The resulting ideals might have new 
linear equations, which we can use to eliminate some of the variables. Geometrically, 
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it means that the projective closure of the motion cm-ve Ca C SE{3) C P^^ spans 
only a subspace P" of smaller dimension. 

The easiest available further information are the Betti numbers of the free reso- 
lution of the coordinate rings of the Ca C P". These numerical invariants determine 
the degree and the arithmetic genus of the curve. Frequently one can deduce the 
connectedness of the curve, which implies irreducibilty in case of smoothness. For 
further information about Betti numbers, we refer to Eisenbud [EOS] and Schreyer 
[S86] . 

The majority of curves of degree 6 leads to curves C C P^ C P^^ whose syzygies 
in P-'' have Betti table 

1 _ _ _ _ 

- 9 16 9 - 

- - - - 1 

In particular, the motion curve C is an arithmetically Cohen-Macaulay curve in 
these cases. The Cohcn-Macaulay property of C implies that if C is smooth, then 
C is also irreducible. In summary, C is an elliptic normal curve of degree 6 in these 
cases. Unfortunately, none of these curves is the reduction of a curve defined over a 
real number field; see Geifi [2008] for an explanation and a unirational description 
of this family. This shows that in general one cannot deduce the existence of real 
exceptional mechanisms from the existence of corresponding families over finite 
fields. 

In case of degree 12 we obtain curves C C P^ with Betti table 

1 - - - - - 

- 10 16 9 - - 

- - 9 16 10 - 

----- 1 

Again the motion curve is arithmetically Cohen-Macaulay, hence irreducible if 
smooth. In this case C is a curve of genus g = 7 in its canonical embedding. 
The fact that these curves occur in their canonical embedding is the main reason 
why we did not work in the Study quadric SE{3) C Q C P^ instead of the naive 
embedding. 

By Schreyer [S86j , the special shape of the Betti table implies that either C has 
a plane model of degree 6 with three (possibly infinitesimally near) double points, 
or that C is bi-elliptic. Indeed, the projection onto S0{3) C P^ yields a plane 
model: 

The projection of C to 50(3) C P^ yields a curve C of degree 12 which is the 
intersection of a Veronese surface with a cubic hypersurface F^: 

C" = n c P^ c P^ 

In particular, C has a plane model C C P^ of degree 6. A general plane sextic 
curve of geometric genus 7 has two different plane models of degree 6. The second 
one is obtained via a Cremona transformation through the 3 non coUinear double 
points. In our case however, the three singular points lie on a line L = {£ = 0}, 
and we have a distinguished model. By adjunction theory, the rational map 

C C C P^ 

to the canonical model of such a sextic is defined by the 7-dimensional space of 
cubics through the double points. Since the double points lie on a line, this space 
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is spanned by 6 quadrics multiplied by i and a further cubic, which intersects the 
line in the three points. We will use all of this information in our construction. 

3. Construction procedure 

Consider Euler parameters on 50(3), i.e., the parametrization by unit quater- 
nions, or, more conveniently, the identification of = iS'0(3) C via the 2-uple 
embedding 

P^-^P^ (so,si,S2,S3) ^ (^(s),t) eP^ 

defined by 

hl-sl-sl + sl -S0S1 + S2S3 -S0S2-S1S3 

A{s) = -sosi - S2S3 -si + sl- s'^ + sj -S1S2 + S0S3 
\ -S0S2 + S1S3 -S1S2-S0S3 -sl-sl + sl + . 

and t = sl + sl + sl + s|. 

In particular, we note that every curve in SO{3) C P^ has even degree. This 
explains why even degree motion curves are more likely than odd degree curves. 
Odd degree can only occur if C intersects the 2-dimensional vertex within the 6- 
dimensional cone SE{3). 



Utilizing the special geometry of the genus 7 curves in SE{3) as a hint, we will 

construct exceptional mechanisms in three steps. We say that a curve C C SE{3) 
supports the leg defined by p'^^^ , q^^^ if the equation eg*^*) is satisfied for all points 

(A, b) e C. 

(1) Choose a plane C through the point (0:0:0: 1) correspond- 
ing to the identity in SO{3) and choose three of the desired six legs 

, , . . . , , qi3) freely. 

(2) Compute the family of plane sextics C" C P^ C = 50(3) together with 
an extension to a birational map C ^ C C SE{'3) C P'^^ such that the 
equations eg'^^eg'^' and eq'^^^ are satisfied for all points on C. It turns 
out that this family depends unirationally on three parameters. 

(3) Choose the parameters such that C supports further three legs. 

To start, we can write the equation eg^*) homogeneously in terms of Euler 

parameters as 

2i2(pW,gW) + 2(ApW,f9W) -2(^p('),fe) + 2(fe,ig«)- II b f==0. 



Since C C 5£^(3)nP^ should be the canonical embedding, the entries of A, t and the 
hi are given by cubic forms on our P^ where the aij and t have the common factor 
i, and the equations evaluate to an equation of the plane sextic. The difference 
eg'*-* — eg*^-'-' has £ as common factor. Division by £ leaves us with a system of 
equations of degree five on the plane. We may regard this as a linear equation 

/bo\ 

[ao ai a2 0) ^ 

for 60) ^ij ^2 and t where the coefficients ai and (3 are quadrics and a quartic which 
depend on the choice of the two legs. Taking two differences, we obtain a system 
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of equations 



ao ct2 P 

a'o a[ a'2 13' 



fbo\ 

hi 

V) 



and any solution substituted in one of the original equations will give a plane sextic 
which supports the first three legs. 

Step 1. We will choose the first three legs and a plane V{s) = c through the 

point (0:0:0:1) corresponding to id <E 5*0(3) with integral coordinates. Then 
S = Q[si, S2, S3] = Q[so, • ■ • ) S3]/ (s), the coordinate ring of P^, is defined over Q. 



Step 2. The solutions {hoMM^^) G S'^ 
kernel of the homomorphism © 5* (—2) 



(B S{—2) are elements of degree 3 in the 
^ 52(2) defined by 

"1 1^' J 

For general choices (at least experimentally), the kernel is the image of a map 
5(-2) © S'(-3) © S{-A) ^ 53 © S{-2) 



a' 



ai 
a' 



defined by a syzygy matrix 



( ao 60 co\ 

ai b'l ci 

a2 b'2 C2 

whose entries are homogeneous forms ai,6-,Ci of degree 2, 3 and 4 and a linear 
form £. In particular, L = V{i) C F"^ is uniquely determined by the choices in 
step 1. For b we have a three parameter family of solutions bi = b^ + gai where 
g = xsi + ys2 + zss is a linear form with three yet unknown coefficients x, y, z. The 
family of sextics is now given by 



/ = 2t2(p(i),g(i)) + 2{Ap^^\tq'^^'>) - 2{Ap^^\b) + 2(6, tq^^)) - bl ~ bj 
which is a homogeneous polynomial in si, S2, S3 with coefficients in Q[x, y, z] 



bl 



Step 3a. We now turn things around and ask for the condition that the curve 
C C SE{3) defined by / supports three further legs. Equations for a leg (p, q) are 



obtained by asking that / and fp 



,Q,x,y,z][si,S2,S3] with 



fpg = 2t^{p,q) + 2{Ap,tq) - 2{Ap,b) + 2{b,tq)- \\ b f 

are proportional as homogeneous polynomials in si, S2, S3. 

This yields an ideal / C Q\p,q,x,y,z] = Q\po,pi,P2,qo,qi,q2,x,y,z] whose 
zero set is the collection of legs. Saturating with the ideal of the known legs, gives 
a simplified ideal J C Q[p, q, x, y, z]. A Grobncr basis computation shows that the 
zero loci X = V{J) C is a surface of degree 18. We are interested in values 
X, y, z such that 

X n V{x — X, y — y. z — z) 

consists of three points. In other words, we are looking for a triple point of the 
projection 

X Z (Z A^, (p, q, .X, z) ^ {x, y, z). 
Since X is a surface, we expect a finite number of triple points. 
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Step 3b. The ideal J is generated by linear and quadratic polynomials in p, q with 
coefficients in Q[x,y, z]. We take these polynomials and all multiples with any of 
the variables p, q. The coefficient matrix with respect to monomials in p, q up to 
degree 3 gives us a matrix with entries in Q[x, y, z]. This matrix drops rank by 3 
at the triple point loci of the projection, we utilize this to compute the triple point 
loci. 



4. Existence over a number field 

For three legs with integral coordinates of moderate height, the computation 
of the ideal J can be done over the integers. But to compute the triple points, that 
is the ideal of minors in Step 3b, is out of reach over Q. Instead we pass ffist to 
various finite fields Fp where this computation is possible. The degree of the triple 
point locus is 210 for nearly all primes. 

It turns out that among the coefficients x, y, z the coefficient z plays a particular 
role: Scanning through the number of triple points on the planes defined over Fp, 
that is z — z for z S Fp, we find for most primes four distinguished values z G ¥p 
which contain a lot of triple points. Three planes contain 28 triple points with 
coordinates in Fp, the fourth 14 triple points. There are 112 further triple points. 

It turns out that X is always reducible: Its projective closure X is 

X = p2 u u u y c P^. 

The residual surface F is a singular conic bundle over P^. One can check this 
computationally using the special shape of the syzygies of Y. We do not explain 
this computation here. Each of the image of a contains 28 triple points, while 
the plane corresponding to the 14 triple points arises as follows: 

Two of the conies get identified under the projection, and the projection of Y 
has 14 triple points along the image conic. We decided to go for a point on this 
conic. The other family of triple points will most likely lead to number fields of 
degree 112, and a further extension of degree 6 for the coordinates of the three 
preimage points. We do not study this second family any further here. 

Using the LLL-algorithm [LLL82j or Wang's rational conversion [W81| . we 
can find the equation in Q[x,y,z] for the special plane, for the conic within this 
plane, and finally for the septic, which cuts the conic in the 14 triple points. Thus we 
obtain number fields Q C K C L of degrees [K : Q] < 14 and [L : K] < 2, where K 
contains the coordinates x,y and the possibly quadratic extension L Z) K contains 
all coordinates of the legs. The degree [L : K] < 2, because the preimage point of 
the triple point, which does not lie on the union of the two conies is distinguished. 
The conies themselves are just conjugated for most choices of three initial legs. 

The specific example below is obtained by choosing the legs for i — 

1, 2, 3 as 




and by taking P^ — V{so — si) C P'^ for the plane. We have choosen this example 
out of many others since the plot of the plane sextic, looked nice to us, see Figure 
4 below. 
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From our computations we obtain the equations for the coefRcients 

_ 90 
^- if' 

the conic 

2 214 377 2 266242 499174 45419195 

O — X H xii v H X v 

^ 181 181 25521 25521 ^ 1199487 

and the septic, which we display below. 

S = 40529281858604847902260745656467021xy^ 
-13003502555830472808716754914961492t/'' 
+1823658874752524908287468326069696796x2/5 
-1345103945703671198677310366396663022y« 
+21057817082494339849424039333426043801a;y* 
-21097516746760090326807604348970652732i/'5 
+104146801533327515594541432822537162912x1/3 
-147061992679082503452275851430863761654y'' 
+261417542057521780350077326987208938863x1/2 
-5735738433410027305745258115342224982362/3 
+137627172885093280917715345045703910468xy 
-11741274978159976338253495611531691168261/2 

-1021528567727081989434308426376376661X 
-15806295004170301422736264965300029948681/ 
- 1752162935565295197876926071913693437570 . 



The septic equation describes rational function of degree 7 of y and 

substituting this into the quadric, yields the polynomial which defines the field 
extension K. 

5. Real solutions 

To prove the existence of a non-empty open set of initial choices as in Step 1 
such that the number fields K C L have a real embedding, it suffices to establish 
this for a single example. Figure 3 shows a plot of the quadric and the septic of our 
specific example from Section |4l There are 6 real intersection points, out of which 
five are visible in the plot. We choose 

{x,y) « (2.2072318327235894538,-7.4995700484259023117). 

The corresponding plane sextic has indeed three double points, all real with real 
tangents, as we can see in the plot of the curve in an affine chart given in Figure 4. 
We conclude from the picture that this specific mechanism will have two disjoint 
connected components over R. 

Next, we compute the remaining legs [p^*-*, q^*-*] for i = 4, 5, 6 using the ideal J 
from Step 3a. They turn out to be real as well. The approximate values are 

- /3.1691819518341l\ /4.06043490691685\ " 

2.10535216460007 , 2.08171150266153 
. \0.09679237835995/ \3.29773800483749/ . 
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Figure 3. A quadric and a septic with 6 real intersection points 
(5 are visible). 




Figure 4. Plane sextic with three double points on a line 

- /2.21062671304078A /3.10187966812352\ " 

1.14679692580674 , 1.1231562638682 
. \~2.1941220884757/ \1.00682353800184/ . 

and 

- / 3.5973027547094 \ /5.9479638679496l\ " 

3.5973027547094 , 4.40897636315916 . 
. Vl. 79865137735471/ V 1.4739819339748 / . 
However, in this position the mechanism is rigid, because the point corresponding 
to id G SE(3) does not lie on the curve of motion. To obtain a start positions 
on each connected component, we apply the motion given by one element on each 
connected component. Of course, the q^*' do not change. The coordinates of the 
p*-*-* for a start position on left hand side component of in Figure [5] are given by the 
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Figure 5. An exceptional Stewart-Gough platform. The upper 
left is the rigid identity position in the assembly space. The re- 
maining pictures show the mechanism in different positions on both 
connected components of the motion curve. 



6-tuple of points with aproximate coordinates: 



0.155967547996383 \ 

= I -0.330589644955676 
-0.915719119302685/ 



^-1.89049475348576\ 

-2.26381911010301 , 
-1.95235593308989/ 



3.02807585569042 \ 
p(^) = ( -1.02097882111392 
,-1.43950977817185/ 



3.54839453008461 \ 
0.988949245364636 , 
,0.195967470045645/ 



3.34246862886166 \ 
p(^) = ( -1.05120251881891 
-1.50147789095895/ 



3.83683978692307 \ 
^(6) ^ I 3.14661620544245 . 
^0.948829119787298/ 
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The right hand side component has a start position: 



/ 0.651863003594838 \ 
= -0.291286328884768 , p^^^ = 
\ -1.84428994846972 / 




/-0.0897065626926926\ 
= 0.326427753808606 , p(^) = 
\ 0.99622125383145 / 



/-0.125381784762398\ 
= 0.392226097496023 , p'^^^ = 
\ 1.30926149281988 / 



Using a singular value decomposition version of path tracking due to Charles 
Wampler, one can trace out the motion curve numerically. 



The family of mechanisms constructed depends on the 3 • 6 parameters of the 
first three legs, and 2 parameters for the C = 50(3) which passes through 
the point id € SO{3). Hence altogether we have 20 parameters, and the corre- 
sponding family M C A^^ has codimension 16. Considering equivalence classes up 
to translation, rotation, and dilation subtracts 7 dimensions, thus we have a family 
of dimension 13. Coming back to the experimental exploration, we did not expect 
to find points on component of codimension much larger than 12. 

Problem 1. Explain why we were so lucky to find examples over F3 neverthe- 
less. 

One possible explanation is that there arc indeed many different components of 
degree 12 mechanisms in A^^. Each could contribute with the same amount to the 
number of Fp-rational mechanisms. For example, the symmetry among the six legs 
is broken. The two legs lying on the conies, seem to have always the same length. 
Already in case of two distinguished legs, we obtain (2) = 15 different components 
in the space of 6 ordered legs. However, this is still not enough to account for the 
difference to the data from the experimental exploration. 

The generators of the syzygies of the matrix 



Problem 2. What is the reason for the low degree syzygies? 
The surface X = U U U y of extra legs has a very special structure. 

Problem 3. Explain, dim X = 2 and the special structure of the surface X. 

In particular, explain the conic bundle structure of Y. 

Problem 4. Why are two conies of Y identified under the projection X 



6. Open questions 



\a'Q a[ a'2 (3' J 
in Step 2 are of lower degree than the syzygies of a general map 




S^(SS{-2) 5^(2). 



Zc A3? 



Problem 5. The legs with coordinates on the two conies have the same length. 
Why? 
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It is not clear whether or not the rcsuhing family of mechanisms is unirational. 

Problem 6. Determine whether this family of mechanism is unirational. 

For a non-unirational irreducible variety M defined over Q or M the set of real 
points A/(M) can have several connected components. 

Problem 7. If the family of mechanisms is not unirational, determine the 
number of its connected components over the reals. 

Since we do not have a unirational parametrization of the family, it is not easy 
to change the mechanisms continuously. Perhaps, it is possible to set-up a path 
tracking approach following the ideas of Sommese and Wampler [SW05j . 

Problem 8. Vary the parameters of the mechanisms using homotopy methods! 

Problem 9. Compute the other 112 triple points and a corresponding mech- 
anism! 

Problem 10. Study other families of mechanisms! 
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